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We introduce and analyze a novel mean-field model for polariton condensates with velocity depen¬ 
dence of the effective polariton mass due the photon and exciton components. The effective mass 
depends on the in-plane wave vector k, which at the inflection point of the lower polariton energy 
branch becomes infinite and above negative. The polariton condensate modes of the new mean-field 
theory are now sensitive to mass variations and for certain points of the energy dispersion the po¬ 
lariton condensate mode represents fractional quantum mechanics. The impact of the generalized 
kinetic energy term is elucidated by numerical studies in ID and 2D showing significant differences 
for large velocities. Analytical expressions for plane wave solutions as well as a linear waves analysis 
show the significance of this new model. 


Introduction.- About two decades ago the fractional 
Schrodinger equation (FSE) was discovered as a mathe¬ 
matical extension within the Feynman path integral for¬ 
malism by transposing Brownian with Levy-type paths 
mm- This generalization of the fundamental equation 
of single-body quantum mechanics has given rise to new 
intriguing mathematical structures and it forms the base 
of fractional quantum mechanics [THS]. The FSE incor¬ 
porates the concept of an intrinsically nonlocal fractional 
kinetic energy, 

(-A)V(r) = .F-i(|k|2*^(/)) 

= / |kPV'(k)e*'^-'dk, (1) 

while the linear SE is the special case 5 = 1. J-{f) = 
/(k) = /(r) dr denotes the Eourier transform 

of /(r) = ^ 2 ^ /(^) ^ On the other hand the 

concept of velocity dependent mass is well established in 
solid state physics [7] suggesting a possible route for the 
implementation of fractional quantum mechanics or even 
more complex kinetic energies as will be shown in this 
letter utilizing polariton condensates. 

To introduce the concept of generalized kinetic energy 
we turn to the solid state system of polariton Bose- 
Einstein (BE) condensates - macroscopically occupied 
single mode states that highlight properties of fundamen¬ 
tal quantum mechanics ranging from quantum harmonic 
oscillators [HE] to interference [iniiii] while providing 
control over key system parameters mna. We show 
that the type of kinetic energy in Schrodinger like models 
is of fundamental importance for the modes and partic¬ 
ularly for non-equilibrium polariton condensate behavior 
at different locations of the dispersion. Polariton con¬ 
densates have kinetic energies of the mathematical form 
of a Eourier-multiplier, jF~^{g{'k)P'{f)), ^(k) is a real¬ 


valued function associated with the two branches of the 
polariton dispersion through [7] , 

5(k) = ^^--= 9{k) = edlE^,vik) (2) 

with k = |k| for k G and d = 1,2. The two en¬ 
ergy branches of the dispersion curve Fl,u(^) vary signif¬ 
icantly over k and the kinetic energy (§ depends on the 
k of the injected or spontaneously populated condensate 
polaritons generally in a non-parabolic way. In fact lo¬ 
cally fractional kinetic energies can be implemented due 
to the velocity dependent mass m(k) that modifies the 
parabolic dispersion accordingly - e.g. a modification of 
the polariton condensate wave function due to effectively 
negative mass was recently shown experimentally m- 

In this letter the whole spectrum of the lower polari¬ 
ton branch is considered while taking the dynamical be¬ 
haviour into account. We clarify the role of the gener¬ 
alized kinetic energy as it is particularly important for 
implementations above the inflection point and because 
several mathematically different forms of the kinetic en¬ 
ergy have been used in similar scenarios while neglecting 
the inherent mathematical inconsistencies of the corre¬ 
sponding predictions as secondary effects [I7HI9]. Cur¬ 
rent models catch aspects of the condensate wave func¬ 
tion at localized k but the concept introduced here in¬ 
corporates the mean field treatment for more extended 
wave packets in k space while being the more accurate 
description even for localized wave packets. Concepts 
such as energy relaxation can be included in the new 
PDE [TTl [19]. Numerically we find that a time split¬ 
ting Eourier pseudospectral method Eoiin] can be used 
to generate converging solutions, a method that will be 
presented in more detail in a later work. 

Theoretical Background. - Polaritons are quasiparticles 
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consisting of excitons and cavity photons within semicon¬ 
ductor micro-cavities which obey Bose-Einstein statis¬ 
tics El and thus the potential to condense into a single 
particle mode [22]. Excitons are coupled pairs of elec¬ 
trons and holes of oppositely charged spin-half particles 
in a semiconductor held together by the Coulomb force 
between them m- Excitons interact with light fields 
[23] and can form integer spin polariton quasiparticles 
in the strong coupling regime that are confined to the 
micro-cavity [24|. As polaritons are 10^ times lighter 
than rubidium atoms [22], condensation is observed in 
CdTe/CdMgTe/GaAs micro-cavities [M] |22j [25] and re¬ 
cently even at room temperature in flexible polymer- 
based structures [26l |27]. The basic Hamiltonian tak¬ 
ing the interaction between the cavity light modes and 
excitons into account is stated in El US]- By diagonaliz¬ 
ing this operator one gets the lower and upper polariton 
eigenvalues El, 

^L,u(k) = i(^Sc(k) + S.(k)T 

T^iE,{k)-E,{k)y+4h^niy (3) 

The dispersion of the cavity photon is Ec(k) = 

where kx,c denotes the 
orthogonal part of the 3d wavevector and using the 
notation ujq = 27rcN/{nL)^ the effective cavity pho¬ 
ton mass rUc = 27rhnNj {cL) and the exciton energy 
Ex{p) = ^ + P^/(2ma,), which can be assumed as con¬ 
stant close the centre of the polariton dispersion. Eor 
our investigation, we set Ex = 1.557 eV, in accordance 
with recent results presented in [25]. n is the refraction 
index between the cavity mirrors, c the speed of light, 
L is the cavity spacer length and N the number of the 
quantized z-mode orthogonal to the k-plane. The polari¬ 
ton mass of each branch depends on the effective exciton 
mass rux ^ 0.1 — Irrie and the mass of the cavity photon 
rric ^ 10~^me with the electron mass rUe and is pro¬ 
portional to the inverse of the second derivative of the 
dispersion®. 

In Eig. [^(a) we show the two dispersions of the polari¬ 
ton, the exciton and cavity photon and in (b) the kinetic 
energy (prefactors) for constant, fractional and velocity 
dependent mass. Eig. [^(b) shows that locally fractional 
and generalized quantum kinetic energies are present due 
to the varying curvature of the effective mass - an exam¬ 
ple is given for s = 5/6, which approximates the bottom 
of the polariton dispersion at /c ^ 0 to a higher accu¬ 
racy than the parabolic dispersion. The effective mass 
switches sign from positive at k < to negative at 
k > fenf. In-between it becomes infinite on a circle cen¬ 
tered around the origin at /c = 0 in the 2D k-plane - the 
inflection point k = k[^f ^ 1.3952/im“^. While models 
of coupled PDEs separating the photonic and the exci- 
tonic fraction have been discussed previously [211120] , we 


present now a unifying approach for the mean-field of 
condensed polaritons resulting in a single non-local PDE 
as realistic model of generalized fractional quantum me¬ 
chanics in a highly controllable solid state system. 
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FIG. 1. (a) Energy dispersions of the upper and lower po¬ 
lariton energy branch (solid lines correspondingly), the cavity 
(parabolic) and exciton (almost constant) dispersions (dashed 
lines), (b) Kinetic energy prefactors g{k) defined in Eq. [^ 
including velocity dependent mass (solid line) compared with 
constant mass (dashed line) and fractional kinetic energy k^^^ 
(dotted line). The sign of the kinetic energy switches for ve¬ 
locity dependent mass. 
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FIG. 2. (Pseudocolor plot) (a) \5'ilj-k\ including velocity de¬ 
pendence of the effective mass for parameters r’cc = 0 and 
Vy = 1 in k-space. (b) \S'ip-k \ for the same parameters but for 
constant mass. From blue - zero magnitude to red stronger 
magnitude. 



Generalized state equation.- Phenomenologically the 
condensate wave function is governed by a complex 
Ginzburg-Landau-type partial differential equation El 
IMISIIITIIII], which includes the effects of polariton self¬ 
interactions, polariton-reservoir interactions and non¬ 
equilibrium properties such as gain and decay of conden¬ 
sate polaritons. An accurate quantum theory of polari¬ 
tons is provided in m- While in the mean-field regime 
the spin of polaritons can become apparent through cir¬ 
cular polarization of the driving light source or TM-TE 
splitting and even spontaneously El, we assume the spin 
coherent case for the introduction of the velocity depen¬ 
dent mass concept. We define the velocity dependent 
mass of the lower branch © (see 0) and include it math¬ 
ematically setting q{r) = E~^[g{k)]. Thus the kinetic 

energy becomes ' /) = (Hi) ' E{f)) = 
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up to a constant and with ^ denoting a convolution. 
Consequently the new polariton state equation resem¬ 
bling a coherent driving scheme [Til [5^ reads as follows: 


idttp{r,t) = {1 - irj) ■ h- q*tp{r,t) + iP{r,t) - 

+ [a\'ip\'^+ V{r,t)+u]ip{r,t) (4) 

The frequency of the bottom of the lower brach is denoted 
cj, (a is the self-interaction strength, V (r, t) is an external 
potential, while 7 is the loss rate of polaritons due to their 
decay. The coherent pumping field is [Il[32] 

P(r,t) = Po(r)e^^^-^e-^^^' (5) 


with Po denoting the pump profile amplitude, ki denot¬ 
ing the 2 d pump wave vector and cji the pump frequency. 
Fig. [^(b) shows that locally close to the particular k un¬ 
der consideration fractional kinetic energies emerge and 
thus the presented model includes as a special case 
a feasible implementation of (driven) fractional quantum 
mechanics. Via ki we can choose experimentally at which 
particular k the condensate wave function is formed on 
the dispersion in k-space. Alternatively with incoherent 
driving schemes one could control k of the condensate 
by adjusting the spot size, which determines the final 
velocity of condensed polaritons [ 12 ]. 

Plane waves of the veloeity dependent model. - Let us 
now present analytical plane wave expressions respecting 
the velocity dependent mass of the polariton system. We 
consider a Id scenario, assume homogeneous pumping, 
no external potential V = 0 and set as ansatz for the 
stationary solution ' 0 (x,t) = 


= iPo- 


2 "T' 

uji--uj fjo{x) 


(6) 


As discussed in the supplemental material [28] by set¬ 
ting — /^i) = 'ipoS{k — K[) we obtain three analytical 
solutions [28]. The physical solution is 




+ _ 
0 — 


(1 T T i6) 

“ 2^t^{-27ia^Po + ^/4{-3aa - Mbaf - 729a4p2)i/3 
(1 - iy/3){-27ia^PQ + ^4(-3aa - 3ihaf - 729a4p2)V3 

6 • 2V3 


using the abbreviations a = uj\ — uj - ^and b = 

7 / 2 . The density tends to zero for Pq -> 

0 corresponding to no pumping of polaritons into the 
condensate and it increases monotonically with Pq. The 
solutions are modified by the velocity dependent mass (§ 
and the constant mass case is obtained by substituting 
niLfk) nfiL. For the plane wave scenario respecting 


m{k) implies including its value at ni of the dispersion 
(see Fig. (b)), which modifies the magnitude of the 
wave function (or the luminosity of the micro-cavity). In 
particularly the solutions respect the negative mass effect 
and note that the Id plane wave solutions can be trivially 
extended to 2d. Nonlocal effects can be expected in more 
general pumping schemes as shown below. 

Linear waves analysis for the generalized driven 
system.- The linearized elementary excitation equation 
around the stationary states without trapping and under 
coherent pump with k[ = 0, for '0(r, t) = y^exp(—icjit) 
is stated in oin] and the eigenvalues of the Bogoliubov 
operators for the plane wave modes respecting the m{k) 
dependence are the generalized Bologliubov dispersions 
of excitations as stated in the supplemental material [28] . 
The signs of the dispersion correspond to the positive and 
negative Bogoliubov branch and as expected the presence 
of the velocity dependent mass adapts the excitation en¬ 
ergy by a non-parabolic k dependence of the effective 
kinetic energy. 

The driven Gross-Pitaevskii equation without ex¬ 
ternal potential, energy relaxation (with effects discussed 
mathematically in [33]), absorbing 00 into the phase and 
setting h = 1 reads 

i—= qi^ — ijfjiP. (8) 

To consider the movement of linear waves with velocity 
V we make an ansatz of the form = (fo Sfj, where 
00 represents the unperturbed part solving the general¬ 
ized driven GPE and a small perturbation on top of 
it. We can think of 0o as one of the solutions derived 
in the previous section. By inserting this ansatz in 
introducing a chemical potential /i (to be identified later) 
and dropping terms of order we get the Bogoliubov 
equation for the perturbation as discussed in more detail 
in the supplemental material [28] . We denote the Fourier 
transform of the pump as P = i f e~^^'^(P — 70o) set 
w.l.o.g. 00 = \/n exp(ido) and choose the constant phase 
of the pump such that 0 o ^ \/n and define /i = an. 
Thus in the reference frame moving with the linear wave 
with velocity v on top of the condensate 0 o we get the 
explicit solution 

^ P/X + P*(i7 + k ■ V - g(|k|) - /j) 

72 + 2/X5r(|k|)+g(|k|)2-2i7k-v-(k-v)2‘ 

(9) 

(T)This solution is a natural extension of the equilibrium 
atomic BEG solutions presented in [36l [37] and those for 
constant mass discussed in [T4l|30]. In Eig. |^we compare 
the solutions due to constant and velocity depen¬ 

dent mass. Eor the sake of simplicity we assume P = I, 
7 = I, /r = I and g(k) is given as described above. The 
results indicate a significant difference in the linear wave 
condensate dynamics which particularly will be investi¬ 
gated numerically in more detail in the following section. 
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FIG. 3. (a) Red line corresponds to the M including the m{k) 
effect over time t with [t] = ps. The blue line corresponds to 
the model with constant mass. Parameters are outlined in 
[3^ . (b) The red line corresponds to M including the m{k) 
effect. The blue line corresponds to the model with constant 
mass m- Units are [ki] = fim ^ while the M corresponds to 
the total occupation number. 


Numerical experiments.-NJe consider for Eq. [^the 2d 
scenario of harmonic trapping V (r) = kr‘^ with strength 
k and a Gaussian pump profile Po{x) = ^dexp — , 

where A is the amplitude, d denotes the position of the 
pump and a its width. In Fig. (a) we show a com¬ 
parison of the total mass of a perturbed condensate wave 
function ||'0||2 = J ^2 \ip\‘^dT = M (L^-norm) as it varies 
in time. To generate the perturbation we start with an 
initial condition 'ipi{x^y) = exp(—and evolve it 
according to @ for m = const, and m = m(k). While 
constant mass implies oscillations of the -norm/total 
density at a given instance of time the velocity depen¬ 
dence of the mass acts as a damping term as shown in 
Fig. I (a). For larger times both models induce con¬ 
vergence to a stationary state of different L^-norm. To 
elucidate the differences in the total mass of the conden¬ 
sate wave function we show in Fig. (b) a comparison 
for different ki As ki increases M monotoni- 
cally increases for the m{k) model while it decreases for 
the classic model, hence offering an experimentally fea¬ 
sible test of the new theory Furthermore in Fig. 

I^we show a comparison of the 2d density distributions 
of stationary states for different ki = v^/2(a,a) with 
a G {0,10.38}, [k] = and accordingly chosen uj[ 

in a mexican hat potential (i.e. harmonic trapping with 
a tight gaussian at the center ED- We observe that 
the density cloud/luminosity for velocity dependent mass 
contracts as k increases (see as well [28]) - a behavior 
analog to attractive atomic BEG in traps [38] . Instead of 
negative/attractive self-interactions, the negative mass 
induces a relative sign between the kinetic energy and 
the still repulsive self-interactions leading to the ob¬ 
served contraction consistent with the experimental re¬ 
sults in m, The diameter of the ring shaped condensate 
increases with k for the m{k) model Fig. [^(a), (b) while 
it decreases for the classic theory (c), (d). In addition 
we simulate 0 reduced to Id and obtain the qualitative 


results presented [28] again showing the increase of total 
mass. Furthermore we provide in [28] more detailed nu¬ 
merical results for the case with and without harmonic 
trapping. 



FIG. 4. (Pseudocolor plots) Gomparison between the 2D wave 
functions governed by the m{k) respecting model in the upper 
line - (a) density and (aa) phase for a = 0 and (b) density 
and (bb) phase for a = 10.38. Gorrespondingly in the lower 
line, density (c), (d) and phase (cc), (dd) stemming from the 
classic model. Blue corresponds to lower densities and red 
is associated with higher densities. Parameters chosen as in 

m- 


Conclusions.- We have identified the polariton con¬ 
densate wave functions with those of fractional quantum 
mechanics by respecting the velocity dependent mass of 
polaritons in the governing PDE. More generally because 
the k dependent kinetic energy of the polariton conden¬ 
sate deviates significantly from the parabolic form new 
phenomena could be observed. Remarkably for A: ^ 0 
a fractional nonlinear Schrodinger-type equation is the 
more accurate model compared to the classic parabolic 
nonlinear SE type models previously used. Via a fea¬ 
sible coherent pumping scheme - driving the polariton 
to condense at a chosen single point of its dispersion - 
one can effectively switch between different points of the 
energy dispersion enabling to test the effects of velocity 
dependent mass. The importance of the in-plane momen¬ 
tum for the emerging polariton condensate was shown 
in the explicit analytical expressions stated and further 
the linear waves analysis suggests significant different dy¬ 
namical behaviour. Consequently numerical simulations 
for Id and 2d scenarios for which a novel time splitting 
Fourier pseudospectral method has been developed re¬ 
vealed evident differences in the predictions for the po¬ 
lariton condensate to the classic results. The dynamics 
show a suppression of total density oscillations due to the 
velocity dependence of the mass, the total mass increases 
for larger k while classic mean field models predict a loss. 
The latter is a phenomenon is a feasible test of the theory 
presented here. In addition condensates forming above 
the infiection point show attractive-type density profiles 
in accordance with the observations in [16] . While the co¬ 
herent driving scheme utilized in this paper pre-defines 
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the phase and suppresses the spontaneous emergence of 
excitations such as vortices or dark and bright solitons in¬ 
coherent driving schemes may reveal interesting pattern 
formation in future works. 
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